Creation of Stable Multipartite Entangled States in Spin Chains with Defects 
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We show how defects in a spin chain described by the XXZ model may be used to generate 
entangled states, such as Bell and W states, and how to maintain them with high fidelity. In the 
presence of several excitations, we also discuss how the anisotropy of the system may be combined 
with defects to effectively assist in the creation of the desired states. 
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I. INTRODUCTION 

Applying the superposition principle to product states 
of composite systems leads to one of the most striking 
properties of quantum mechanics: entanglement. What 
used to be a subject within the foundations of quantum 
mechanics is now a resource for the development of new 
technologies, playing a central role in quantum compu- 
tation and quantum communication. 

The generation of highly entangled states is a prereq- 
uisite for quantum teleportation, various quantum cryp- 
tographic protocols, and is also needed in quantum com- 
putation. In the bipartite setting, maximal entanglement 
appears in the Bell [or also called Einstein-Podolsky- 
Rosen (EPR)] state (1/V2)(|10) + |01>) Q. Diir et al 
showed that there are two different kinds of genuine 
tripartite pure state entanglement: the maximally en- 
tangled Greenberger-Horne-Zeilinger (GHZ) state [3j and 
the so called W state. The later is the state of three 
qubits that retains a maximal amount of bipartite en- 
tanglement when any one of the three qubits is traced 
out. It is written as (l/v / 3)(|100) + |010) + |001)). Here, 
we show how Bell and W states may be created in spin 
chains with defects. 

Spin chains provide ideal systems for the study of en- 
tanglement. They are naturally used to model quantum 
computers: the two states of a spin-1/2 particle corre- 
spond to the two levels of a qubit and the exchange in- 
teraction corresponds to the qubit-qubit interaction. We 
focus here on a spin chain described by the XXZ model, 
where two different kinds of interaction are identified: the 
Ising part, proportional to the anisotropy coupling, and 
the XY part, which can be used to create entanglement 
between two or more qubits 0, . 

In principle, the qubit level spacings may be control- 
lable j(| , allowing the creation of entanglement between 
precisely selected qubits. as it turns out, the anisotropy 
of the system can also be used to our advantage, for states 
with different numbers of neighboring excited qubits are 
not coupled and a large system may then be treated as 
consisting of several uncoupled small chains |5|, |jj . 

The paper is organized as follows. The model is de- 



scribed in the next section. Section 3 discusses how de- 
fects can be used to create Bell and W states when the 
system contains just one excited qubit. Once the de- 
sired state is generated, we outline a way to preserve it 
with high fidelity for a long time, which requires quickly 
detuning the defects involved in the process. Section 4 
shows how the anisotropy may be used as an extra tool 
toward the creation of entanglement. 



II. THE MODEL 

A spin chain with nearest neighbor coupling is consid- 
ered. The Hamiltonian describing the system has two 
parts: one corresponds to the Zeeman energy of each 
spin, Hd, and the other is related to the spin-spin inter- 
action and is given by the XX Z model, Hxxz- In units 
where % = 1, we have 
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Above <7 z ' + '~ are Pauli matrices. There are L sites and 
we deal with a periodic (or closed) chain, that is, sites 
n + L and n are the same. Each site n is subjected to a 
magnetic field in the z direction, giving the energy split- 
ting e n . In terms of qubits, e„ is the level spacing of 
qubit n. A spin pointing up corresponds to an excited 
qubit, or an excitation. The parameter J is the hopping 
integral and A is a dimensionless parameter related to 
the anisotropy coupling. The non-diagonal term -Shop 
is responsible for propagating the excitations and delo- 
calizing the states. The diagonal term a n a n+1 gives the 
Ising interaction. It is only relevant when at least two 
excitations are present, being therefore associated with 
many-body problems. We set J and A > 0. 
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We note that some recent works have discussed entan- 
glement in anisotropic systems with impurities 0, 0, , 
but their anisotropy differs from the one in Hamilto- 
nian The hopping part of the Hamiltonian con- 

sidered here can be equivalently written as ifhop °c 
JxCn a n+i + Jy a n a n+i- Here J x = J y = J, but in the 
models cited above, the degree of anisotropy comes from 
the difference between J x and J y . The anisotropy in our 
case originates from the extra Ising interaction. In terms 
of entanglement, few studies have been developed with 
this model 0, , though this is the model that best 
describes some proposals of quantum computers [|| . 

In a homogeneous chain, all qubits have the same 
level spacing e, while a disordered system has at least 
one site with a different energy e n = e + d n , which is 
called a defect. In principle, complete control of the 
qubit level spacings is available, which allows the 'cre- 
ation of defects'. We set d n > and assume that e 
largely exceeds J, A and d n , so the ground state of the 
system corresponds to all spins pointing down. In what 
follows, we count energy off the ground state energy 
So = — X)n=i £ «/2 + LJA/A, i.e., we replace in Eq. (Q) 
H — ► H — Eq . 

The states \aia2--.oir,) , where = or 1, correspond 
to the computational basis. To address the different 
states of the system we use a notation that is common 
when determining the energy spectrum of spin chains 
with the Bethe ansatz ^jj. A state with a single ex- 
citation on site n, that is | |i| 2 •■■ in-itn-U+i •■■ II), or 
equivalently |0i0 2 ...0 n _il n 0„+i...0£), is simply written 
as 4>{n). A state with two excitations, one on site n and 
the other on site m, is <p(n, m), which is a simplified no- 
tation for | J.xJ.2 ■■• Tnln+i ••• Tm ••■ II), or equivalently 
|0i0 2 ...l„0„+i...l m ...0 L ). 

In the system considered here, the z component of the 

total spin, S z — Yln= l "-'«> ^ s conserve d> hence states with 
different number of excitations are not coupled and the 
Hamiltonian effectively describes uncoupled subsystems. 
The diagonalization of a certain block corresponding to 
N excitations leads to L\/[N\(L — N)l] eigenstates, such 
that the kth eigenstate is written as 



L 

^ a (fc) (ni,n 2 , ...,n N )\<p(m,n2, ...,n N )). 

n 1 <ri2<...<n N = l 

ILL SINGLE EXCITATION 

Suppose that there is only one excitation in the chain 
and assume that all qubits have the same level spacing e, 
except for some defects that have all the same level spac- 
ing e + d. By choosing d much larger than the hopping 
integral J, we guarantee that states with one excitation 
on a defect do not couple with states that have the ex- 
citation on the other sites. In terms of energy spectra, 



we form two well separated bands. It is then possible to 
treat the spin chain as two smaller and decoupled chains. 



A. Two defects: Bell states 

Let us first examine the case where there are simply 
two defects with d^> J, placed on sites n\ and n 2 . The 
chain is broken into two, one with L — 2 qubits, whose 
state energies lie within the band £\ ± J, where £\ = e — 
J A, and the other with just two qubits, whose eigenstates 
are Bell states and have energies ~ £\ + d. An excitation 
created on one defect will hop between the two defects 
with a frequency determined by their distance. To find 
this value, we need to compute the energies E± of the two 
eigenstates i/j± = 4j[^>(m) ± 0(w 2 )] of the defect chain. 

If the two defects are nearest neighbors, n\ and ri\ + 1, 
they are coupled to first order of perturbation theory. We 
diagonalize a 2x2 submatrix whose diagonal elements are 
£i+d and off-diagonal elements give the effective hopping 
integral J e s = J/2, so E± = E\ + d ± J/2. If the defects 
are next-nearest neighbors, ri\ and ri\ + 2, they are cou- 
pled to second order of perturbation theory. Now the di- 
agonal elements are £i+d+ J 2 / (2d) and the off-diagonal 
elements are J c s — J 2 /(4c?), leading to states with much 
closer energies E± = £ 1 + d + J 2 /(2d) ± J 2 /(Ad). The 
larger the separation between the defects, the closer will 
be the energies of the two eigenstates. 

An excitation initially prepared on site n\ will have 
probability 

l + cos[(E + ~£-)t] 

^(ni)W = 2 ' > 

to be later found on the same site. The oscillation period 
of the excitation between the defects is inversely propor- 
tional to the energy difference of the two eigenstates. It 
therefore depends on the number fi of sites between the 
defects as T M = T (2d/J)' J , where T = 2n/J. 

Bell states are created when PM ni \ — PM n2 ) = 1/2, 
that is, at every instant ts = nk/[2(E + — J5_)], where k 
is an odd number. The more separated are the defects 
the longer is the wait for a maximally entangled state to 
be created, but in principle we can entangle even remote 
qubits. 



1. Detuning procedure 

Once the desired state has been generated, the two 
defects should be detuned for the state to be maintained 
unchanged. A scheme for a possible detuning procedure 
is shown in Fig. ^ The level spacing of one qubit is kept 
constant, while the other increases according to a certain 
function of time S(t), which depends on the experimental 
setting and on the fidelity we want to attain. At the 
end of the process, the energy difference between the two 
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defects should be much larger than their effective hopping 
integral. 



1 
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FIG. 1: Scheme for the detuning procedure for two defects on 
sites rii and 712 . 

The detuning should start at t = ts and, apart from 
constants, it is well represented by the Hamiltonian 
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dctu 
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Jeff 



(3) 



In Fig. [3 we show the results of a linear detuning, 
S(t) — Dt, and a quadratic one, S(t) = Dt 2 . Two defects 
with d = 10J coupled to second order of perturbation 
theory are considered. The larger the ratio D/J e g, the 
closer the generated state will be to the Bell state at later 
times. Contrary to what happens when creating states, 
now more separated defects make it easier to preserve the 
state with high fidelity. 
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concurrence Q and global entanglement 0, 0, Il7| . 
can be affected by defects. They both vary from (un- 
cntangled state) to 1 (entangled state), but concurrence 
is a measure of bipartite entanglement, while the other 
is able to quantify multipartite entanglement. 

It is straightforward to conceive a system where all 
eigenvectors are Bell states and therefore have maximum 
individual (as well as the average over all states) con- 
currence. An example is a chain with an even number 
of qubits L, where pairs of sites have the same level 
spacings, but differ from the others by more than the 
coupling strength between them. In this situation, all 
eigenvectors are Bell states each involving a particular 
pair of resonant sites. In terms of multipartite entan- 
glement, designing a distribution of defects that would 
lead to large global entanglement is in general more 
demanding. The definition of global entanglement is 
2 — (2/L) J2n=i ^fiL where p n is the reduced density 
matrix for qubit n [l5l.ll6f. For the example described, it 
is maximum only for a two-qubit chain, decreasing with 
the chain size as 2 — 2/L — 2(L — 2)/L. Moreover, systems 
described by the XX Z model may not be ideal, because 
states that are known to have large global entanglement, 
such as the GHZ or products of Bell states, cannot be 
created. This is a consequence of the conservation of the 
total spin in the z direction, which prevents the coupling 
between basis states that have a different number of ex- 
citations. 



B. Three defects: W states 

To create a W state with one excitation, three defects 
are necessary. The procedure is very similar to the one 
described before. Let us assume three nearest neighbor 
defects, m, n% = ni + 1, and = m + 2, with d~^> J. 
The three eigenvectors of the defect chain are a good 
approximation to the actual corresponding vectors ob- 
tained with the total Hamiltonian. They can be simply 
computed by diagonalizing the submatrix 



£1 + d J/2 
J/2 £ x + d J/2 
J/2 £1 + d 



FIG. 2: Time evolution of the probability to find the excita- 
tion on site m during the detuning of two defects placed on 
ni and m + 2. The left panel has S(t) = Dt and the right 
panel has S(t) = Dt 2 . The detuning is initiated at the mo- 
ment a Bell state is created. The effective hopping integral 
is J e ff = J 2 /{Ad), d — 10J, and the values of D are shown in 
the figure. 



2. Measures of entanglement 



giving 

^ a = |[0(n x ) + V2(j)(n 2 ) + 4>(n 3 )} E a =£ 1 + d+ J/y/2, 
1* = 775 l-Hni ) + <P(n 3 )} E b 



ipe 



V2 



£i+d, 

(ni) - V2cb(n 2 ) + c/>(n 3 )} E c =£ x +d- J/y/2. 



The W state is generated by preparing an initial state 
with an excitation in the middle site ri2, since this defect 
is coupled to the other two to first order of perturbation 
theory. The probability for the excitation to be found on 
this same site later in time is 



Several measures of entanglement have been intro- 
duced. Here we briefly comment on how two of them, 



)(*) 



1 + cos[(£ Q - E c )t] 



(4) 
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while the probabilities for sites n\ and are the same 



1 - cos[(£ a - E c )t] 



(5) 



For the W state we need P^ ni )(t) = P^ n2 )(t) 
P(f>(n 3 )(t) = 1/3, which happens at every instant 



(-l) fc arccos(-l/3) + 27r(fc- |fc/2]) 
tw = p p , (6) 

where fe is an integer. 

A possible detuning procedure for the defects for t > 
tw is the following: the level spacing of defect nz is kept 
fixed, while the level spacings of rii and na increase in 
time according to the functions 8i(t) and 6%^), respec- 
tively. The Hamiltonian describing the process is given 
by ' 
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For the three defects to become out of resonance, we can 
use for instance, Si(t) — D\t, 5 2 {t) = D 2 t, and D\ ^ 
D 2 . The faster the detuning is performed the higher the 
fidelity for the desired state. Moreover, the final result 
becomes better if D 2 > D\ instead of D\ > D 2 . An 
example is shown in Fig|3 




FIG. 3: Time evolution of the probability to find the exci- 
tation on site n\ and on site m + 1 during the detuning of 
three neighbor defects. The probability for the third defect 
placed on iii + 2 is very close to the probability for m . The 
detuning is initiated at the moment a W state is created. For 
the left panel D\ = 10J and D2 = 100J, for the right panel, 
Di = 50J and D 2 = 500J. 



IV. MULTIPLE EXCITATIONS 

When more than one excitation is present in the chain, 
the Ising part of the Hamiltonian starts having an effect. 
If the anisotropy parameter A is large, i.e. J A 3> J, 
several well separated energy bands appear, each corre- 
sponding to states with a different number of excitations 



next to each other. The case of two excitations, for in- 
stance, is schematically shown in Fig.0] The states where 
the excitations are far from each other have energy in- 
side the band 2£i ± 2 J, while the states where they are 
next to each other have a much higher energy and form 
a much narrower band 2£i + JA + J/(2A) ± J/(2A). 
Such bound pairs are coupled to second order of pertur- 
bation theory, taking much longer to propagate through 
the chain if compared to the free excitations. The two 
bands are not coupled and can be treated separately. 
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FIG. 4: Energy spectrum of the considered spin chain with 
two excitations and one defect with JA >d> J. 

In addition, suppose that a single defect ni is added to 
the chain with two excitations, such that JA ^> d ^> J. 
Two new bands appear [see Fig. 0] . One band has states 
with the two excitations far from each other, but one 
trapped on the defect, they have energy 2£i + d ± 2 J. 
The other band is made of just two (Bell) states having 
the largest energies. Here, the two excitations are next 
to each other and one of them is on the defect site. We 
have, 



i>± = ~ ni ) ^ ^{ni,ni + 1)] 

E± = 2£i +d + JA + — 



(8) 
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J 



4A 4(JA + d) 4(JA + d)' 

Advantage may be taken of the large anisotropy to 
create entanglement. Preparing an initial bound pair on 
the sites {n\ — l,n\), for example, will lead to Bell states 
at every instant tBP = 2( JA + d)[7r/2 + kn]/ J 2 , where k 
is odd. As before, detuning site n\ — 1 or site ni+1 should 
preserve the maximally entangled state. Moreover, since 
states (n\ — 1, n\) and (m, n\ + 1) are coupled to second 
order, the interaction strength is small and the detuning 
rate required to keep the created Bell state with high 
fidelity does not need to be very large. Several other 
possibilities for Bell and W states where use is made of 
the defects and anisotropy have been discussed elsewhere 
110 
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V. CONCLUSION 

We have investigated how bipartite and tripartite 
states can be generated in strongly anisotropic spin-1/2 
chains described by the XXZ model by controlling the 
Zeeman energy of the spins [or equivalently, by manip- 
ulating the qubit level spacings]. Even though the in- 
teraction between the qubits is always on, it is effective 
only between states with very close energy. Selected reso- 
nant defects can then be used in the creation of entangled 
states. Once generated, the maintenance of the desired 
states depends on how fast the defects involved in the 
process can be moved away from resonance. 

Our interest in the XXZ model emerges from its rele- 
vance to some proposals for quantum computers @, but 
we now intend to extend the method developed here to 
other models, especially those that they may be used to 



create multipartite states with large global entanglement. 
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